Effect of Landau Level Mixing on Braiding Statistics 
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We examine the effect of Landau ievei mixing on the braiding statistics of quasiparticles of abelian 
and nonabelian quantum Hall states. While path dependent geometric phases can perturb the 
abelian part of the statistics, we find that the nonabelian properties remain unchanged to an accuracy 
that is exponentially small in the distance between quasiparticles. 



In 2+1 dimensions, quasiparticles may exist that obey 
nontrivial braiding statistics [lj. A well-known exam- 
ple occurs for the v = 1/m Laughlin quantum Hall 
state, where taking one quasiparticle adiabatically clock- 
wise around another accumulates a Berry's phase e 27 "/ m 
which is independent of the details of the path, but de- 
pends only on its topology (i.e., that it went around the 
other particle) . Since the braiding statistics is comprised 
only of a phase, this case is known as "abelian". In 
more complicated "nonabelian" quantum Hall states [2], 
the ground state manifold can be expressed as a vector 
space, and taking one particle around some others applies 
a unitary matrix to this vector space, with the unitary 
matrix depending only on the topology of the path and 
not on its geometric details. This independence from 
details can in principle allow highly robust quantum in- 
formation processing and is fundamental to the concept 
of topological quantum computation 

It it known, however, that certain effects can introduce 
corrections to the braiding statistics of abelian quan- 
tum Hall states @, H Most seriously, when Landau 
level (LL) mixing is taken into account, the long range 
Coulomb interaction can introduce corrections that are 
power law in the distance between quasiparticles. If sim- 
ilar corrections were to occur for the nonabelian braiding 
statistics of nonabelian quantum Hall states, it would 
present a fundamental barrier to the idea of topological 
quantum computation The purpose of this paper is 
to show that this does not occur. To be precise, we as- 
sume that we are given a quantum Hall wavefunction, 
fully in the lowest Landau level (LLL), which is in a well 
defined topological phase (with all its associated topo- 
logical properties) and we then examine the effects of LL 
mixing. We find, for LL mixing not too large, to accuracy 
exponentially small in the distance between quasiparti- 
cles, only the abelian phases are perturbed — which then 
does not present a problem for the idea of topological 
quantum computation [lj]. 

At first glance, it might seem natural that nonabelian 
statistics would be better protected from perturbations 
than abelian statistics, as it is well known[6] that non- 
abelian statistics can only occur in a discrete set of vari- 
eties (which cannot be perturbed continuously) whereas 
there is a continuum of possible abelian statistics. How- 
ever, this so-called Ocneanu rigidity 6] of nonabelian 



statistics is not sufficient to show that the nonabelian 
statistics are preserved under perturbation, as it possible 
that the result of taking a quasiparticle around a loop 
may end up being non-topological (i.e., may depend on 
the details of the path and not just its topology), and 
hence the perturbed system may not have a well defined 
braiding statistics at all. Indeed, this is precisely what 
happens in the abelian caseQ where the correction (from 
long range interactions and LL mixing) to the phase ac- 
cumulated by taking one quasiparticle around another is 
dependent on the distance between the quasiparticles. 

We start by reviewing how one calculates geometric 
phases and topological statistics associated with trans- 
porting quasiparticles in the LLL [3, 0). For nonabelian 
quantum Hall states, in the presence of quasiparticles, 
the ground state is described by an arbitrary superposi- 
tion of several wavefunctions which we call "blocks" 
(for the abelian case we have only a single block). If 
there are N electron coordinates z,; and M quasiparticle 
coordinates w a , we write these blocks as 
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where £ is the magnetic length, e* is the quasiparticle 
charge, and we choose to work in holomorphic gauge 
where <j> is holomorphic in all of its arguments. Although 
4> is single valued in the z coordinates, it may be multi- 
ple valued in the w coordinates (due to branch cuts) , and 
this multiple valuedness is why there are multiple blocks 
4>i. We assume that these blocks are ortho normalized 
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If wavefunctions are generated as conformal blocks of a 
conformal field theory (CFT) with an appropriately cho- 
sen background charge, then the form of Eq. [1] in holo- 
morphic gauge appears naturally [H, 0|. Orthonormality 
(Eq. [2|) for conformal blocks is then often assumed[l|, 0] 
so long as quasiparticles are well separated, although this 
has only been proven for a few special cases (HQ. 

For nonabelian quantum Hall states, the effect of adi- 
abatically transporting a particle w around a closed loop 
is to apply some unitary matrix to the space of blocks 
— > Uij This matrix is the product of the explicit 
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monodromy matrix M. (which gives the explicit change 
in the wavefunction due to branch cuts when moving the 
coordinate w) and a Berry's matrix B given 



B = V exp 
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where r indexes the position of the quasiparticle, and V is 
the path ordering operator llj . Writing the quasiparticle 
position in holomorphic and antiholomorphic coordinates 
we have d/dr = (dw/dr)d w + (dw/dT)dw- Noting that 
4>i is holomorphic, applied to acts only on the 
gaussian factors to give (^i\dw\^j) — —8ij{e*/e)(w/At 2 ). 
Similarly, using Eq. [3] we have d w (*S? i\& j) = so we 
can replace (vPil^l^j) by — (3^(^1)1^) so that d w acts 
only on the gaussians. We then easily obtain the abelian 



Berry's matrix B — Sije n where 

1 ={A/e){e*/e) 
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where A is the area surrounded by the path. The fact 
that this is proportional to the area enclosed means that 
it is a geometric phase, and is therefore not considered 
part of the braiding statistics. Thus, the braiding statis- 
tics come entirely from the explicit monodromy of the 
wavefunction in holomorphic gauge. (For wavefunctions 
generated by CFT, the nonabelian statistics is thus the 
monodromy of the conformal blocks, assuming Eq. [2]). 

The point of this paper is to extend this calculation to 
the case where LL mixing is included as a perturbation. 
Our starting point is to assume that the unperturbed 
state (fully in the LLL) is indeed a topological state of 
matter with all of its associated properties [lj. In partic- 
ular we will assume (a) that if quasiparticles remain well 
separated then (to accuracy exponentially small in the 
separation between quasiparticles for large enough sepa- 
rations) blocks cannot be distinguished from each other 
by any topologically local measurement [l2j| . In addition, 
we note (b) that since the system is in a gapped phase, 
any sufficiently small change in the LLL Hamiltonian will 
not alter the topological phase of the system in the LLL, 
and in particular will not destroy property (a). This 
is consistent with the fact that for many quantum Hall 
states there is a broad range of LLL Hamiltonians that 
produce essentially the exact same wavefunctions. 

Let us now be a bit more precise about assumption (a). 
Given an operator O that acts nontrivially on electrons 
only in a region 1Z, we say that the operator is "topologi- 
cally local" if the region 1Z consists of disconnected (well 
separated) regions each of which surrounds only a sin- 
gle quasiparticle or no quasiparticles. Assumption (a) is 
equivalent to the statement that for any topologically lo- 
cal O we must have (^i|0|^j-) = CS^ for some C (which 
may be a function of the w coordinates). This assump- 
tion is crucial to the idea of using quantum Hall states 
for topological quantum computation, as the inability of 
any local operator to distinguish blocks is precisely the 



property that makes the space of blocks highly immune 
to decoherence. We emphasize again that this statement 
has not been proven in general, but is frequently assumed 
to be true. Note that if O is any sum of topologically lo- 
cal operators, then (^j|0|^) = CSij remains true, and 
we think of this as also being topologically local. 

To understand the effects of LL mixing, we first de- 
fine the interaction operator V = J dqu(|q|) p(q)p(— q), 
where p(q) is the density operator at wavevector q and 
v(q) = 27re 2 cxp(— q£ v )/(eq) is the Coulomb interaction 
(with e the dielectric constant) which we have cut off 
at some finite short length scale £ v to avoid singulari- 
ties. (This cutoff may be provided in a physical sys- 
tem by finite well width, for example). It is conve- 
nient to work with a coherent state basis [HI for sin- 
gle particle wavefunctions which we write as ^R,n(r) = 
|BL, n) where n is the LL index and R is the guiding 
center coordinate. Although these coherent states are 
non-orthogonal, we still have the completeness relation 
1 = (2Tre 2 )- 1 J2nI dR|R,n)(R,n|. We define annihi- 
lation operators for an electron in such a coherent state 
cr,« = J dr i J 9R ! „(r)^(r) and using completeness obtain 
■0(r) = (2n£ 2 y 1 J2 n f dR^ R ,„(r)c R , n , where tp is the 
usual electron annihilation operator. The density opera- 
tor can then be written as 
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A key observation here is that the matrix element is local: 
|(Rn|e !q ' r |R'?i')| reaches a maximum of order unity or 
less at R — R' = z x q£ 2 (where z points normal to the 
plane) and decays exponentially for all |R — R'| > q£ 2 . 
Since we have cut off the interaction at |q| ~ l/l v this 
means that in a single application of the interaction V ~ 
p(q)p(— q), no electron guiding center coordinate can be 
moved more than a distance of order I 2 /l v . Note that 
even though V is long ranged, it is a sum (or integral) of 
topologically local operators. 

We now define an operator Po to be projection of all 
electrons to the LLL, and define V — V — PqVPq to 
be the interaction, with the LLL part removed (i.e., V 
is the part of the interaction that we want to treat as 
a perturbation). The unperturbed Hamiltonian will be 
Po = K + PqVPq with K the kinetic energy, and the 
full Hamiltonian is then H = K + V = Ho + V. Given 
an eigenstate \9) of the full Hamiltonian, we define the 
LLL piece, l^ ) = Pol*)- It turns out to be convenient 
to work with a normalization where (5 ,0 |5' ) = 1. We 
thus write \^! u ) = \^°) + l^ -1 ), where the superscript u 
means "unnormalized" and we construct the normalized 
wavefunction |*) = Z x l 2 \^ u ) with Z" 1 = (* u |^ u ). 

We write the Schroedinger equation as |\P U ) = (E — 
Po) -1 ^!*") and then project with 1 — P to obtain 
I*- 1 ) = &(E)V(\V°) + I*- 1 )) where we have defined 
G{E) = (1 - P )(P - Po)" 1 = Ek \k)(k\/(E - e k ) 
where the sum is over only states k (with unperturbed 
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energy £&) which are not fully in the LLL. Note that 
the energy denominator is (at least) of order the cy- 
clotron energy hio c . We then obtain the full wavefunction 
|*") = X(E)\V°) where X(E) = ET=o(G(E)V) k . We 
have thus written a perturbation expansion for the full 
wavefunction \^ u ) in terms of its LLL component |\I/ ). 

We now find an effective theory for wavefunctions 
fully in the LLL. We can rewrite the Schroedinger equa- 
tion now as HX\^°) = EX] 1 ® ) which, projected to 
the LLL, can then be recast as H eff \^ Q ) = E\^°), 
with the effective Hamiltonian H e ff = -^o + SHq where 
SHq = P VXP = P [VGV + VGVGV + . . .]P . Here 
SHq is a LLL interaction that includes all of the effect 
of integrating out higher LL's exactly. The perturbed 
energy is thus given by E = Eq + (^°\SHq\^ ) where 
co = ( x f ,0 |i/o| x f' ) is the unperturbed energy. 

Our initial assumption is that the unperturbed sys- 
tem (H in the LLL) is in a gapped topological phase. 
As mentioned above in (b) since the state is assumed 
to be gapped, perturbing the Hamiltonian slightly with 
the LLL interaction SHq cannot change the topological 
phase of matter within the LLL. Thus, we assume that 
the system defined by Hq + SHq fully in the LLL is also in 
this topological phase, and in particular satisfies assump- 
tion (a) above. We denote the blocks of this perturbed 
LLL systems as \^). From these blocks, we can use our 
above perturbation theory to generate the corresponding 
full wavefunctions \ as described above. What we will 
show below is that the nonabelian part of the braiding 
statistics of the blocks I'F;) of the full multi-LL system 
must be the same as that of the LLL blocks |\&°). 

First, let us re-examine Eq. [5] for the blocks of the 
full system. We define [X^X]„ to be the operator X'X 
expanded to n th order in perturbation theory. (Each or- 
der in perturbation theory is a factor of GV and is or- 
der E cou i om b/hui c — e 2 /(e£huj c ) ). As emphasized above, 
each application of V can only move an electron a dis- 
tance of order £ 2 /£ v (and application of G does not move 
anything at all). Thus, the operator [A^A]„ can only 
move particles a maximum distance of n£ 2 /£ v , and is 
therefore topologically local so long as quasiparticles are 
a distance d > n£ 2 /£ v apart. Since the LLL wave- 
function is in a topological phase, assumption (a) 
tells us that (^°\[X^ X} n \^°) = CS^ with accuracy in- 
creasing exponentially with d for d > n£ 2 /£ v . So long 
as our perturbation theory converges [lj], then for suffi- 
ciently large n the difference between the matrix element 
of [XtX] n and X^X will decrease exponentially quickly 
with increasing n. Thus, for sufficiently large d, we have 
{^\X^ X\^) — CSij with accuracy increasing exponen- 
tially with d. Thus, to this exponential precision, Eq.[2]is 
obeyed by the full wavefunctions \^>i). Furthermore the 
normalization constant (\&"|\I r ") = Z^ 1 is independent 
of the block i. Similarly, using locality of SHq, we show 
that the energy E is also independent of the block. 



We are now prepared to calculate the braiding statis- 
tics of the system with LL mixing. As above, we first 
calculate the Berry's matrix Eq. [3] Here we use the 
wavefunctions = Z 1 ' 2 ^) = Z^X^), We 

must now be cautious because X and Z depend not 
only on the quasiparticle position, but also on the 
energy E. We thus write d/dr — (dw/dr)d w \ E + 
{dw / dr)dyj\E + {dE / dT)dE\ w is where we use derivative 
notation where parameters listed after the | in the sub- 
scripts are held constant. Fortunately, as shown in the 
appendix, i\d E\ tu .w\^ j) = 0) so we can drop this term 
from d/dr. We now proceed as above in the unperturbed 
LLL case. Application of d^E to \^>j) now acts only 
on the gaussian factors of and on Z. We obtain 
= 5 l3 [-(e*/e)(w/U 2 ) + (1/2)9^ log Z]. 
As above we use orthonormality (Eq. [2]) to move 3 w \e to 
act on (^il giving a similar result. We again end up with 
an abelian Berry's matrix with phase 7 + ^7 where 7 is as 
given above in Eq. |4]and #7 = — Im § dw 8 w \e log[Z -1 ]. 
Crucially, there is no nonabelian contribution from the 
explicit Berry's matrix, although the abelian contribu- 
tion has a correction. Thus the nonabelian braiding 
statistics are again entirely given by the explicit mon- 
odromy. However, since both the normalization Z 1 / 2 and 
the energy E are independent of the block, and (by lin- 
earity) the monodromy of X\^) is obviously the same as 
that of 1 4'°), we find that the monodromy of the full wave- 
function is the same as that of its LLL component 
I \E f °) . This is the main technical result of this paper: the 
nonabelian part of the braiding statistics is independent 
of projection to the LLL (subject to the same condition 
that the quasiparticles are far enough apart and correc- 
tions are exponentially small in this distance). Note that 
our perturbative analysis relies on having a cutoff for the 
Coulomb interaction at some finite distance £ v (which is 
indeed true in experimental systems). One can ask if the 
result remains true without this cutoff. While our deriva- 
tion cannot handle this case rigorously, it seems unlikely 
such details of the interaction would be important. 

We note in passing that we can use similar techniques 
to study perturbations of systems completely in the LLL. 
For example, given a solvable Hamiltonian Hq with a 
topological ground state in the LLL (such as the spe- 
cial interaction that gives the Read-Rezayi wavefunction 
as its exact ground state[l5|]), and given a perturbation 
interaction V in the LLL, we can use an analogous argu- 
ment to study Hq + V where the perturbation parameter 
is now V/A with A the excitation gap of Hq. Indeed, 
assuming (a) holds for blocks of Hq, we can show that 
the blocks of Hq + V also satisfy (a) so long as V is appro- 
priately topologically local. This can be used as a rough 
"proof of assumption (b) above. 

We now return examine the abelian phase 7 + ^7 
calculated above. Since Z" 1 = (^°\X^X\^°) = 
1 + (*°|V r G 2 V r |\E' ) + the correction 67 is order 
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[E coulomb/ fuuc} 2 ■ This result apparently contradicts prior 
work [i], Q which suggest corrections at first order. How- 
ever, there is a subtlety that we have not yet consid- 
ered: in the unperturbed wavefunction the param- 
eter w represents the position of the quasiparticle (i.e., 
its complex coordinate is z = w). However, once the 
wavefunction is perturbed, the quasiparticle's complex 
coordinate will generally move to some perturbed posi- 
tion z = w + 5w. Neglecting the second order correction 
6~f, the geometric phase associated with moving the pa- 
rameter w around a given path remains given by 7 in 
Eq. [4] where A is the area enclosed by the path in param- 
eter space. However because the quasiparticle's position 
is perturbed, the physical path of the quasiparticle is not 
the same as the path in parameter space and does not 
enclose an area A. 

We now estimate the magnitude of this effect. Let 
us imagine a quasiparticle at the origin as the source 
of an unscreened long range potential U(r) — ee*/(er). 
We first solve the problem of noninteracting electrons ex- 
posed to this potential. In symmetric gauge, we write the 
single particle wavefunction with angular momentum m 
as ip m = e %mS y(r)r~ x l 2 and we obtain an effective one 
dimensional Schroedinger equation — {fuv c /2)d 2 y/dr 2 + 
u(r)y = Ey where the effective potential is u[r) = 
U(r) + (fiw c /8)[-4m + (4m 2 - l)(£/r) 2 + (r/£) 2 }. Solv- 
ing this perturbed Schroedinger equation mixes LL's of 
the unperturbed basis. The resulting (perturbed) ground 
state wavefunction is very peaked near the minimum of 
the effective potential, which in the long distance and 
weak interaction limit is at position r max w \j2m£ + 5r 
where Sr = I [e* /(2em)][(e 2 /e£)/(fuv c )]. This shift in po- 
sition of each orbital (assuming no other change to the 
multiparticle wavefunction) gives precisely the 1/r 3 cor- 
rections to the density profile of the quasiparticle found 
in Refs. 0, 0|- As a result, any other quasiparticle at 
distance r from the one at the origin will accordingly 
have its position shifted by 5r, although its position will 
still be labeled by the original parameter w. To re- 
turn the quasiparticle to position z — w, the param- 
eter w needs to be changed by roughly —Sr. Taking 
this quasiparticle around the one at the origin then gives 
an extra abelian geometric phase — {2irr(8r) / £ 2 ){e* / e) = 
— (2irl/r)(e*/e) 2 [(e 2 /e£)/(huj c )] in comparison to the case 
without the long range interaction and LL mixing. 

It has also been pointed out that there is yet an- 
other source of correction to the abelian statistics due 
electro-dynamical effects [1, Here, the currents in the 
quantum Hall layer generate a very small amount of ac- 
tual physical flux, and this physical flux binding slightly 
changes both the abelian statistics and the total charge 
of the quasiparticle[3, 5]. However, this mechanism again 
has no way to effect the nonabelian braiding statistics. 



To summarize, we have demonstrated that, to expo- 
nential accuracy with increasing distance between quasi- 
particles, there are no corrections from long range in- 
teractions and LL mixing to the nonabelian statistics of 
quantum Hall states. This result is crucial to the idea 
of using such quantum Hall states for quantum informa- 
tion processing [J. We have also examined perturbations 
to the abelian phase from LL mixing and find several 
contributions, the leading term of which agrees with pre- 
vious works [j, [j| . While these calculations have all been 
done for LLL wavefunctions, it is trivial to generalize to 
the higher LL's, where we suspect that we may have seen 
actual nonabelian states of matter experimentally 

The author acknowledges conversations with N. Bon- 
esteel, V. Gurarie, M. Freedman, J. Slingerland, and 
F.D.M. Haldane. 

Appendix: Throughout this appendix we fix quasipar- 
ticle positions. As in the main text we use locality of 
X*d B X to show (9i\d B \9j) = {9?\Xld B X\&}) = C<%. 
Since = Z l / 2 X\^) is normalized, (^i\d E \^ l ) can 
only be nonzero if it is imaginary. Using d E G — —G 2 we 
show d E X = -XGX + XG so that {^\X^d E X\^) = 
-(^\X^XGX\^) since G|*°) = 0. However, since 
XG = GX* and G is hermitian, we have X'XGX hermi- 
tian. Thus the matrix element is real so i\d E \^ j) = 0. 
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